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Abstract 

We investigate the collapse of a spherical shell of matter in an anti-de 
Sitter space. We search for a holographic description of the collapsing shell in 
the boundary theory. It turns out that in the boundary theory it is possible 
to find information about the radial size of the shell. The shell deforms the 
background spacetime, and the deformed background metric enters into the 
action of a generic bulk field. As a consequence, the correlators of operators 
coupling to the bulk field are modified. By studying the analytic structure 
of the correlators, we find that in the boundary theory there are unstable 
excitations ("shellons") whose masses are multiples of a scale set by the radius 
of the shell. We also comment on the relation between black hole formation 
in the bulk and thermalization in the boundary. 
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1 Introduction 



One of the many exciting aspects of the AdS/CFT correspondence [0, 0, [§] is 
that it contains the seeds for solving the black hole information puzzle in a 
holographic || 0] way. So far, the understanding of black holes in the context of the 
AdS / CFT correspondence is mostly limited to the properties of static solutions. For 
example, one can give a boundary theory interpretation for the Bekenstein-Hawking 
entropy 0, |J , properties of different vacua and their thermodynamic interpretation 
|10| , PI , [13|, [14], However, the static solutions have some rather special 



features, such as the appearance of multiple asymptotic regions*]. To address the 
issue of the fate of information, one would like to study black holes which form from 
collapsing or colliding matterf], and their subsequent evolution to the end. These 
issues are much less clear as discussed in In this paper we take a first step 

towards understanding the collapse of a spherical shell of matter in AdS space. 

The first problem to confront is how to account for the presence of a shell of 
matter in the bulk of an AdS space. In particular, how does one obtain information 
about the radial size of a shell? For a pointlike probe, the scale-radius duality 



21] , [22] , [12] , [23[1 converts its radial distance in the bulk into a transverse scale of its 
image at the boundary. However, in the case of a shell, the radial scale is a bit 
trickier to find because of spherical symmetry. What is the holographic image of a 
shell at the boundary? Since a pure AdS space corresponds to a vacuum state in 
the boundary, and the end result of the shell collapse (if it is a stable black hole 
with positive specific heat) corresponds to a boundary theory at finite temperature 



23, a first guess is that adding a shell to the AdS space corresponds to an 



off-equilibrium configuration in the boundary theory, which then evolves towards 



thermal equilibrium as the shell collapses to a black hole in the bulk (see e.g. fll2" , |25 
for discussion). Further, the scale-radius duality of the AdS/CFT would suggest that 
the initial configuration is somehow peaked at higher energies than the temperature 
of the black hole. In other words, one might hope to associate an initial average 
occupation number spectrum with the initial configuration, with a peak at energies 
corresponding to the initial radius of the shell, and study how the time evolution 
rearranges the spectrum to a thermal one, with the peak moving down in energies 
to stabilize at the Hawking temperature. 

To test the above idea, one would need some way of detecting or computing 
e.g. the average occupation number spectrum in the boundary theory. One would 
then need to study operators which have access to the excited states. For example, 



4 However, see p6| . 

5 Black hole creation in AdS3 has been studied in J!?], [l^, |i~9|j . In these cases the resulting black 



holes are stable. 
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the shell deforms the background spacetime, and the deformation couples to the 
stress-energy tensor of the boundary theory, giving it an expectation value. The 
expectation value gives information only about the total energy, not about the details 
of the distribution. However, one could study higher n-point correlators of the stress 
tensor; these involve more detailed information about the spectrum which one can 
then try to extract out0. 

Another, more direct, way to detect the shell would be to assemble it from matter 
associated with a bulk field coupling to an operator C s h e ii in the boundary. Then, 
the n-point functions of s heii would provide information of the shell itself, without 
a detour to its backreaction to the metric. However, this route obviously depends 
on what the shell is made of and how. 

In this paper, we shall follow a third possible strategy. Since the shell deforms the 
background spacetime, the action for any bulk field in the deformed background is 
modified as well. Since the bulk action is the generating function for the correlation 
functions of the corresponding operators in the boundary theory, the correlation 
functions will reflect the presence of the shell. Even if the expectation values of 
the operators (the one point functions) vanish, the 2-point functions can be differ- 
ent from those evaluated in the vacuum. For example, in the end of the collapse, 
when the boundary theory thermalizes, a two-point function will become periodic 
in imaginary time, even if the operator has a zero expectation value in the thermal 
background. So our strategy is that we shall not try to directly probe the states 
associated to the shell, but we investigate their effect on other operators in the 
boundary theory, and use their correlation functions to obtain information about 
the shell. A related possibility would be to study the properties of a hanging string 



pq , |27| , |28t P9| . In the case of a black hole the boundary theory will be at finite 
temperature, and the force between two quarks will be screened at sufficiently large 
separation. In the picture with the quarks corresponding to endpoints of a hanging 
string, screening occurs when the string breaks and one obtains a configuration with 
a string from each of the quarks going straight down through the horizon of the 
black hole. With a shell the situation is different since the boundary state is not 
thermal and screening is not expected. From the bulk point of view there is no 
horizon and the string between the quarks can therefore not break. Even though 
the distance dependence of the force will be affected by the presence of the shell, 
there is still a long range contribution. 

The simplest two-point function is given by a bulk scalar field <p coupling to an 
operator O in the boundary theory. For example, could be the five dimensional 



6 In general, higher n-point correlation functions are needed to have more complete information 
about a quantum state. For example, in [20| it was discussed how higher n-point functions are 
relevant for understanding the causal structure in AdS/CFT. 
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dilaton in AdSs, and O the operator TrF 2 in the 3+1 dimensional N = 4 SYM 
theory. Also, we will concentrate on the simplest possible situation. We consider 
a spherically symmetric thin shell of matter, and focus on the initial stage of the 
collapse when the shell is moving very slowly. In that case, we can use a quasistatic 
approximation and take the shell radius to be fixed. We solve the field equation of 
4> on both sides of the shell, impose appropriate boundary conditions, and join the 
solutions across the shell. We will then derive the two point function^ 

(0(t,x)0(t',x')) 

for the corresponding boundary operator O. Actually, our calculation will give the 
momentum space Fourier transform of the two-point function, 

G(u, k) = fdt f d d x e- Ml *(O{t,x)O{0,0)) . (1) 



The analytic structure of the propagator gives information on the spectrum of 
excitations of the theory. In the presence of the shell, the two point function will be 
identified as the retarded propagator G rct (uj, k). For a fixed value of k, the retarded 
propagator is analytic in the upper complex u plane. It can be written in an integral 
representation form 

n ( t\ [°° duJ ' A (u',k) 

GretW.fc = / — — , (2) 

where A(u, k) is the spectral density function. It is the imaginary part of the 
retarded propagator, 

A(u>, k) — —2- lmG ret (uj, k) . (3) 

(See e.g. [[Kj for discussion.) The spectral density function can be interpreted as 
the probability^ that an excitation created by the operator O has an energy u and 
a momentum k. Let us now recall some possible analytic features of the retarded 
propagator in a generic QFT. If we try to do an analytic continuation to the lower 
half-plane, and take the limit e — > 0, we may find that the propagator has poles at 
the real to axis. In the boundary theory, these correspond to elementary excitations 
or stable bound states (the latter being possible in an interacting theory). The poles 
of the propagator correspond to delta function peaks of the spectral density function. 
There may also be cuts, signaling multiparticle states^. We may also find poles of 
the propagator in the lower complex plane, corresponding to unstable bound states 



7 Note that in the case of the static shell, the two-point function will depend only on the 
differences (t — t') and (x — x'). 

8 For bosonic excitations, one needs uj > 0. 

9 Of course, this issue played an important part in studies of absorption by branes ]3l| , in the 
path which led to the AdS/CFT correspondence. See for more discussion and references. 
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(resonances). If we denote the location of such a pole by uJi{k) — zTj, the imaginary 
value Tj is the width of the unstable bound state. The closer to the real axis the pole 
is, the longer lived the resonance is. In the spectral density function, the resonances 
appear smeared peaks with finite width. (See e.g. j32|, |3(| for discussion). 

Thus, the analytic structure of the propagator gives information on the stable 
and unstable excitations present in the boundary theory. In general, the operators 
O are composite operators built out of the elementary fields of the boundary theory. 
Hence, we obtain only information about a certain subset of particles propagating 
in the boundary, the stable and unstable composite objects that can be created and 
annihilated by the operator O. 

The key fact is that this gives us enough information to detect the radial size of 
the shell. We discover unstable resonances in the boundary theory, whose masses 
are multiples of a mass scale set by the radial size of the shell. The boundary theory 
is conformal, but it is in an off-equilibrium state characterized by a dimensional 
parameter, so the massive resonances can exist. 

The resonances are discussed in Section 2. In section 3 we discuss the results. 



2 Propagators in the shell background 

We shall now proceed, step by step, to compute the two point function of a bound- 
ary operator in the case that a (quasi) static shell of matter is present in the bulk. 
Although this is an unphysical situation, it can be used as an approximate descrip- 
tion of a slowly moving shell. In 0, 0, it was described how to compute boundary 
correlators of an operator 0(x, t) using classical solutions for a related bulk field 
0(r, x, t) using the AdS/CFT correspondence. The computation method we will use 



is a variant of the one described in the appendix of ||33|| . 

In all cases, the spacetime is asymptotically AdSd+i (d = 2,4), and near the 
boundary the bulk metric reduces to the approximate form 

ds 2 —f^dt 2 + —dr 2 + —dx 2 (4) 

where the boundary is at r — > oo. Let be a scalar field in the bulk, with a mass 
m. In Minkowski signature, the classical equation of motion 

-^d^ (y/gg^duttr, x, t)) - m 2 0(r, S,t) = (5) 

has two kinds of solutions: normalizable and non-normalizable ones. We denote 
the former by < - + - ) (r, x, t) and the latter by (f>(~>(r,x,t); they have an asymptotic 
behavior 

(j) {±) (r, x, t) -> r~ 2h± (r -> oo) (6) 
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where 

h± = ±(d±y/<P + 4m*) = ±±^. (7) 

The starting point is to impose appropriate boundary conditions for the bulk field 
4> at the interior of the AdS space. In the asymptotic region, the field will in general 
become a superposition of the normalizable and non-normalizable solutions, with 
an asymptotic behavior 

<f>{r,x,t) « r' 2h+ (j) + (t,x) + r~ 2/l -0_(t,£) . (8) 

In the above, is the boundary data which acts as a source for an operator 

0(t,x) in the boundary theory. Next, we Fourier transform the variables t,x and 
write the equation in the form 

0(r, u, k) « (r~ 2/l - + r^+GV, k)) 0_ (w, £) , (9) 

where 

= tt4- < 10 > 

<p-{u},k) 

One can check that the two point function in the boundary will be 

(O(ou, k)0(u', k')) = (h+ - h^)5{u + u')5{k + k')G{u, k) . (11) 

Thus, the essential steps of the computation are as follows: one imposes interior 
boundary conditions for the bulk field, writes it as a superposition of the normal- 
izable and non-normalizable modes, then the ratio of the superposition coefficients 
gives the two point function in the boundary. To obtain the answer with a cor- 
rect overall coefficient (h + — hJ), one needs to perform the calculations with more 
subtlety, as discussed in |33|]. However, overall coefficients will not be important 
here, so they will be suppressed in what follows. The field <f>(r, x, t) must also obey 
time boundary conditions related to the fact that one can compute, for example, 
retarded, advanced or Feynmann propagators. In the absence of interactions, and 
with respect to a vacuum, all propagators in momentum space are equal. Otherwise, 
the propagators can be identified by examining their analytic structure. 

As a warm-up, we consider the propagator in the absence of a shell in the bulk. 
The full AdSrf + i manifold is covered by global coordinates, with the metric 

r 2 + B 2 B 2 r 2 

ds = "T* + ^hv dr + h dS1 ^ (12) 

where the radial coordinate r extends from to oo and dQ 2 i _ l is the metric of a 
(d — l)-sphere of radius B. The interior boundary condition is that must be 
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regular at the origin r = 0. The solution is expressed as a linear superposition of 
normalizable and and nonnormalizable modes in [TT| . The propagator ( JIT| ) is simply 
the ratio of the relative coefficients, which are found in equation (37) of Ref. JTT . 
The answer is 

rt v\ - r (~^ r(/l+ + I ^ ± k))T{h+ ± ffcg + k » 

U[UJ > fcj T(u) T{h_ + f ( W + fc))r(/L + f (-a; + fc)) ' 1 J 

where /c = integer /R. The above result is valid when v is not an integer. If v is an 
integer then the propagator can be written in terms of derivatives of T functions. 
In this case it is easier to compute first the spectral function part and then use the 
integral representation @ to reconstruct the full propagator. It is easy to see that 
the propagator fllcf) has poles at uo± k = =p(2n + h+)^ n — 1, 2, . . . corresponding to 
the normalizable modes in AdS space. From the boundary point of view the discrete 
spectrum is possible since the spatial directions are compactified to a sphere S d ~ 2 . 
Accordingly, note that the dispersion relation is not Lorentz invariant. Between 
successive poles, the propagator ( |i~3"D has also zeroes at uo ± k = =F(2n + hJ)^. 

In the limit R — > oo, keeping uo and k fixed, the poles and zeros accumulate and 
they degenerate to cuts in the propagator for uo > k and uo < —k. Further, in the 



limit R — > oo the boundary decompactifies. Applying the limit to (|13|) , we find that 
the propagator reduces to the form 

G(w, k) ~ {uo 2 - k 2 f {uo 2 > k 2 ) . (14) 

Thus, we recover the expected result for a free propagator in an infinite boundary. 
From the bulk point of view one has to take R — > oo with r/R 2 fixed or equivalently 
r — > oo. In this limit the metric reduces to (H). Using this metric the propagator 



(0) follows immediately. The only subtlety is the boundary condition to use in the 
interior since in the large R limit, the point r = corresponds to the coordinate 
horizon of a Poincare patch. Near the horizon any solution can be written as a linear 
combination of an ingoing and an outgoing wave. If uo has a small imaginary part 
it is easily seen (see Appendix) that one has to use an ingoing wave for Im(u;) > 
and an outgoing one for Im(u;) < 0. The use of different boundary conditions for 
lm{uo) > and Im(co>) < produces a cut on the real axis corresponding to switch 
from the retarded to the advanced one. 

Let us now turn to the case of a black hole of radius r + > R in AdS space. Now 
we must impose a boundary condition at the horizon. Again, both solutions of the 
wave equation are regular at the horizon. The general solution is a superposition of 
an outgoing and an ingoing wave: 

4> w A{r — r + Y% + B(r — r + ) -1 ^ {r—>r + ). (15) 
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Imposing the condition that <p be regular at r = r + implies that if Imu; > we take 
the solution with A = and with B = if Imoo < 0. Since the solutions with A = 
and with B = are one complex conjugate of the other, the imaginary part of the 
propagator changes sign when u crosses the real axis. In the case of a BTZ black 
hole of radius r+, the propagator ([[(J) follows fronQ |L3[ [37], |38|, |39|, ||| and is 



r ( u\ v n( _L ' h\ r(l - v) T(h + ~ f(w + k))T(h + - f (a; - fc)) 

G r pt(a;, K) = llllG W + lC, ft = — : r — ; tt- rr — ; 4r~ rr 

reU ' ; e -o v ' ; r(i + i/)r(/i_-f (w + ife))r(ft_-f (w-ib)) 

(16) 

As before, if z/ is integer then the calculation must be done more carefully. For 
example, for v — 1 the result is 

G(u, k) = {uj2 ~ k2) J J (17) 
i e 2T - l e 2T - l 

This propagator can be understood as a thermal propagator of a composite opera- 
tor. Acting on a vacuum, the operator O creates states with right- and leftmoving 
particles of two-momentum k\ = (ui,Ui) and k 2 = (002,-002)- These two-particle 
states have total energy to = lo\ + oo 2 and total momentum k = uoi — uo 2 . With this 
notation the propagator reads 

G(u,k) = ± i J? £ 2 ■ (18) 
J- e t — 1 e t — 1 

Now it is obvious that this is a thermal propagator for a composite operator prop- 
agating two particle states. The Bose occupation numbers arise from the gamma 
functions ||40|| . In fact this propagator follows simply from conformal invariance, so 



it does not give us any information about interactions in the boundary theory, only 
that it is at temperature T ~ r + /R 2 . A similar calculation can be done for the black 



hole in AdS^ using the function obtained in [JL4] . The result is similar to the case of 
AdS 3 although it cannot be expressed in terms of special functions as in fllEl) . 

After this lengthy introduction, we are prepared to compute the propagator in 
the presence of a shell of radius r s . In the interior of the shell, the metric is the AdS 
metric in global coordinates, given by (1121) . The exterior metric is the AdS black 
hole metric, 

where the parameter /1 is related to the ADM mass of the black hole and determines 
the horizon radius r + . The interior boundary condition for the field is the same 



■•For coordinate space expressions including the correct overall coefficient, see also 



8 



as in pure AdS space: regularity at r = 0. The equation (|5p reduces to: 

r l~i d r (r*~7(»W(r, u, k)) + (j^ - J - m 2 ) <j>(r, u, k) = (20) 



with 



1 + r if r < r 



/m^:,:^ if r>r ; m 

and -Is- = 1(1 + d — 2), I E Z> is the eigenvalue of the laplacian on the sphere 
S^ 1 . The solutions of the field equation, <pi for r < r s and <p2 for r > r s, are the 
known solutions in the global and black hole coordinates. The matching conditions 
at r = r s are 

hl=r a = <h\ r=r . (22) 

/iM^iU = /2(r)a r 2 | r=rs (23) 

where the last condition follows from integrating the equation between r s — e and r s +e 
with e — > 0. It also ensures current conservation. Using the matching conditions the 
propagator follows as 

G(tu, k) = 7-t 7-y (24) 

where are the normalizable and non-normalizable modes in the black hole back- 
ground (in three dimensions, they can be found in fll3"|), and <fti is the mode in global 
coordinates that is regular at r = 0. We have performed this calculation, and have 
found a tower of resonances in the boundary theory. However, it will be much 
simpler and more illustrative to discuss a computation which uses some additional 
approximation methods, which are good for a large shell with radius r s ^> r + ^> R. 
In this limit, the metric in the interior reduces to the Poincare metric, so the solu- 
tion in the interior will be a Poincare mode. Correspondingly, the interior boundary 
condition will be mapped to those in a Poincare patch, as we discussed before in 
the context of the large R limit. We will relegate the details in the Appendix, and 
quote just the end result. By studying the propagator in d dimensions, we find that 
it has an infinite number of poles at the complex values 

irr s ( 3 v\ ir s . ( Ann 

u n w — £ n H h- ^ln r - §- . (25) 

n R 2 V 4 2/ 2R? \(d-\)r%) y 1 

As we discussed in the Introduction, these poles are interpreted as unstable excita- 
tions (resonances) in the boundary theory that appear due to the shell background. 
The real part (the mass) is proportional to the radius of the shell. We call the 
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unstable excitations as "shellons". Note also that the width is also of the order of 
the shell radius, hence the shellons have a short lifetime. This is good, since then a 
slowly moving shell is static compared to the lifetime of the shellons, and the results 
are consistent with the quasistatic approximation. 



3 Discussion 



Our calculation shows the existence of unstable excitations, shellons, with properties 
like mass and life time depending on the radius of the shell. Thus, in particular, one 
can detect a difference between spherical shells of different radii but with the same 
mass. It would be very interesting to have a more detailed understanding of how 
the shellons and their properties are to be understood from the boundary point of 
view. 



Another important problem is to let the shell fall freely and follow its evolution 
into a black hole. Let us make some comments on what to expect, but first we need 
a more detailed understanding of the composition of the shell from the boundary 
point of view. We begin with a static shell consisting of a large number of elementary 
particles. A single particle with mass of the order 1/R at rest at radius r s , will 
through holography show up as a blob in the boundary with size R 2 /r s , and total 
energy given by E ~ r s j R 2 . The shell will therefore be a superposition of such blobs, 
and correspond to a state that is out of thermal equilibrium with a homogeneous 
energy density. If we turn to the falling shell, each blob will turn into a bubble 
expanding at the speed of light JTJJ. and a falling dust shell containing several 



particles will then correspond to a superposition of such expanding bubbles. How 
will they interact? 

To understand this, it is useful to consider first a particle falling into a black 



hole. In [|Tj] the shape of a bubble produced by a particle falling into a BTZ black 
hole was given. If one considers the energy momentum tensor as in p5| then the 
shape is given by[^ : 

Too = 2 + 2 ( 26 ) 

(a 2 + (1 + a 2 ) sinh 2 !ii|±£) (a 2 + (1 + a 2 ) sinh 2 ZtSf&j 

where a is a parameter related to the initial position of the particle and r + is the 
radius of the black hole. We see then that in this case the bubble still propagates at 
the speed of light. However, for fixed t, the energy density decays exponentially for 



11 This follows simply by performing the conformal transformation that corresponds to change 
from global to BTZ coordinates in the bulk 
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large x with a characteristic length of 1/T where T ~ r + /R 2 is the temperature of 
the black hole, while if no thermal bath is present, it decays with a power law. If the 
particles consist of supergravity modes, their mass is of order ~ 1/R while that of 
a black hole is ~ N 2 /R and so we need order N 2 of them. In the large N limit the 
gas of bubbles in the boundary will be very dense and we expect that each bubble 
will see a mean field created by the rest. Each bubble will then expand similarly 
as in eq . (p6|) . Alternatively, we can have a dilute gas of bubbles in the boundary 
if the shell consists of k dust particles, each containing of the order N 2 elementary 
particles of mass 1/R. This could be thought of as a shell consisting of k black 
holes (of course in this case it is not a thin shell). In the boundary, one now finds a 
dilute gas of k expanding bubbles. Even when the falling dust approaches r + , it is 
easy to see that the bubbles still remain dilute, unless the boundary equivalent of 
the bulk gravitational interactions are taken into account. The interactions will, as 
the final black hole forms, merge all the bubbles into a homogeneous background. 
Another way to understand this is to note that the holographic scale/radius relation 
is changed when the deviation from a pure AdS metric becomes important, making 
the bubbles even larger. Does this mean that the k small black holes have merged 
into one big black hole? In this context it is important to know the strength of the 
interactions. The interactions between the black holes that make up the shell are not 
suppressed by powers of N, and are therefore part of the supergravity description. 
In units where R = 1, the Newton's constant is G ~ 1/N 2 , but since the mass of 
each of the black holes that make up the shell is M ~ N 2 , [25j the corrections to 
the metric around one of the black holes are GM ~ 1. However, we do not expect 
these collective interactions to result in a final thermal state, since this will require 
interactions between the individual particles that make up the black holes and the 
shell. An interaction between any two elementary particles will be suppressed by 
powers of N, and therefore thermalization is expected to be a 1/N effect. We 
expect this to be a general result, not depending on the detailed composition of the 
collapsing shell. It is intriguing that similar results have been obtained previously in 
studies of thermalization in large N field theories in a different context |42|, |43| . 
In fact the collapse of the falling shell gives interesting information about the non- 
equilibrium processes they have studied. However a precise relation is certainly 
lacking. 



While it is satisfying that we can find information in the boundary theory about 
a collapsing shell in the bulk, it is also obvious that this is just a step towards 
understanding the formation of black holes as a hologram. One of the steps to follow 
would be to go beyond the quasistatic approximation and study what happens as 
the shell is falling rapidly, or what happens when it crosses its Schwarzschild radius. 
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It will also be interesting to study what happens when the resulting black hole is 
small and has a negative specific heat. Another interesting issue would be to add 
charge to the collapsing matter, to make contact with the recently discussed charged 
AdS black hole solutions EL [45]. 
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Appendix 

We are interested in finding the solutions in the region near the radius of the shell, 
where we match the interior solution with the exterior solutions. If the shell is large, 
this region is at r ^> R, where the interior metric takes the approximate form (H). 
Correspondingly, the interior solutions are equal to the Poincare modes, which can 
be expressed using Bessel functions, as discussed e.g. in [fLU] . The situation is now 



analogous to the large R limit which we discussed in section 2. The continuum of 
Poincare modes will correspond to a cut in the propagator. Correspondingly, we have 
a choice of two interior boundary conditions: the solution must look like an ingoing 
or outgoing wave as r — > 0. Which one is appropriate follows by taking the unique 
solution in global coordinates and obtaining its behaviour for large frequencies. 
Alternatively one can reason as follows: an ingoing wave behaves as 



)j(r) « rVe 1 ^, (r -> 0) (27) 



where we set R = 1 as in the rest of the Appendix. If u has a small positive imaginary 
part the function is well behaved for r — > but it blows up if the imaginary part 
is negative. So an ingoing wave is appropriate for the region above the cut, and an 
outgoing wave for the region below. We choose the former condition, so in the end 
we will obtain a retarded propagator. We denote the interior solution by <fii, as in 
section 2. 

Outside the shell the solution to equation fl20|) cannot be expressed in terms of 
known functions (except for the BTZ black hole), so we use a WKB approximation. 
Again we consider the region r r + R. To start with, we rewrite the equation 



42 



(P(]|) in an alternative form. First, we rescale the field as 

4>(r, u, k) = 1 x(r, w, k) (28) 
y/r*-if(r) 

where f(r) = f^ir) = 1 + r 2 — fi/r d ~ 2 . Then, the equation for the function x can be 
written in a form which resembles a Schrodinger equation 

- d„x {±) + V{r)x {±) = (29) 

with a potential 

(d-l)(rf-3) ^./(r) _ 1 /a,/(r)\ 2 1 (d - 1) fl./(r) 
l ' ' " - — * ■ 1 <•< * > h 2 r /(r) 

(30) 



4r 2 2/(r) 4 V /(r) / 2 r /(r) 

^ 2 -(4 + m 2 )/(r) 



/2(r) 

with f(r) as defined above. 

Outside the shell, we want to find the normalizable and non-normalizable solu- 
tions x^i which satisfy the boundary conditions 

Now r ^> r + , and the potential reduces to 

, , v 2 -\ u?-k 2 

= — ^- , (32) 



where v = yd 2 /A + m 2 and k 2 = 1(1 + d — 2) (where Z is a nonnegative integer). In 
this case the equation can be solved in terms of Bessel functions, and we find the 
normalizable and non-normalizable modes to be 



X (±) (r,uj 



, k) » 2±T(1 ± u)(u 2 - k 2 Y v ^J ±v ( VUJ2 k2 ) (33) 



However, in the region (d 2 — l)r 2 /4 <C uj 2 — k 2 , namely for large frequencies, better 
approximate solutions are obtained using the WKB expression 

X = [Aexp(z J y/f^jdr) + B exp(-i J y/V^)dr)^ (34) 

with 

v( r ) = u2 - fcV " 4 m) 
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Matching the WKB functions with the Bessel functions, the following expressions 
are obtained for the normalizable and nonnormalizable modes: 



X {± \r,u,k) « J- \ 2^T(l±u)(u; 2 -k 2 )^ x 
V 7T {V(r)) l / A 



X cos 



7T 7T 



Vuj 2 - k 2 ( r% 2u 2 - k 2 N 

r [ + r d 2(uj 2 -k 2 ) ) T 2 U 4 



(36) 



in the region (d — l)r 2 /4 <C u> 2 — k 2 and r ^> r + . These expressions will be used to 
find the propagator using (pip. 

Thus, for large frequencies, we can use the approximate solutions (^) outside the 
shell. In the large frequency limit, the Poincare mode 0i, discussed above, reduces 
to the form (we can focus on k = 0) 



i 



r—e 1 ' (37) 

up to some constants which cancel in the propagator. Inserting the above expressions 
to the formula (|2~4D , we obtain the propagator (at k = 0) 



r x 2l/ r(l+i/) 2 „ae~^ ( 'e- m 2 + /? e * ( V"^ 
v ; r(l-i/) ae-^ (+) e +Mr 2 +^ l+) e- ,T s 



with 

t(=F) 



c<j / ri \ zr 



i 1 + 73 1W + T ( 39 ) 
r s V (d+ l)r d J 4 v ; 



a = (l^fii-^ (40) 
2 r? r. 



r 



A useful check is that for r + = the last fraction in the propagator is a constant, so 
the propagator correctly reduces to the expected form in pure AdS. The propagator 
has poles when the denominator vanishes. For r s ^> r + it follows that there are 
poles at the values of to which satisfy 

00 3 v i n f-, 4iu; r d \ . 

— = -7T H In 1 + — r §- . (42) 

r s 4 2 2 V (d-l)r s r d + j v ; 

In the leading approximation, we drop the "1" inside the brackets in the last term. 
We obtain 

7tr s ( 3 v\ ir s , / Ann r d \ , . 

^?(" + i + !J-sHm^ ' (43) 
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where n is a (large) positive integer. Notice that the imaginary part is negative since 
r s > r + . 

Finally, we would like to emphasize we have also performed a calculation using 
the exact solutions in global and black hole coordinates, with regularity at r = as 
the interior boundary condition. In that case, the resonances are found by examining 
the imaginary part of the propagator, the spectral function. They are seen in the 
envelope that modulates the amplitudes of the spikes associated with the elementary 
excitations (the poles at the real values of uj. The envelope has periods of hills and 
valleys, and by examining the position and the width of the hill one can deduce 
the position of the complex poles. In this case the calculations were performed 
numerically. In this way one obtains a better approximation for the location of the 
poles, but the qualitative results remain the same. The approximations used in the 
above give accurate enough results to reach the main conclusions in the case of large 
shells. 
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